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Partial  Coordination  Numbers  in  Binary  Metallic  Glasses 


DANIEL  B.  MIRACLE,  KEVIN  LAWS,  OLEG  N.  SENKOV,  and  GARTH  B.  WILKS 

A  critical  analysis  of  measured  partial  coordination  numbers  for  binary  metallic  glasses  as  a 
function  of  composition  shows  a  large  scatter  of  ±1.5  but  clear  trends.  The  current  work  uses 
two  topological  models  to  predict  the  influence  of  relative  atomic  size  and  concentration  on 
partial  coordination  numbers.  The  equations  for  partial  coordination  numbers  derived  from 
these  two  models  can  reproduce  measured  data  within  experimental  scatter,  suggesting  that 
chemical  effects  on  local  structure,  although  present,  may  be  relatively  small.  Insights  gained 
from  these  models  show  that  structural  site-filling  rules  are  different  for  glasses  with  solute 
atoms  that  are  smaller  than  solvent  atoms  and  for  glasses  where  solute  atoms  are  larger  than 
solvent  atoms.  Specifically,  solutes  may  occupy  both  f]  and  y  intercluster  sites  when  the  solute- 
to-solvent  radius  ratio  R  is  less  than  1.26,  but  only  P  sites  can  be  occupied  by  solutes  when 
R  >  1.26.  This  distinction  gives  a  simple  topological  explanation  for  the  observed  preference  for 
binary  metallic  glasses  with  solutes  smaller  than  solvent  atoms.  In  addition  to  structure- specific 
equations,  simplified  phenomenological  equations  for  partial  coordination  numbers  are  given  as 
a  convenience. 
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I.  INTRODUCTION 

The  number  and  type  of  atoms  that  surround  a  given 
central  atom  is  the  simplest  description  of  atomic 
structure  in  metallic  glasses.  This  first-neighbor  consti¬ 
tution  dictates  the  number  and  type  of  atomic  bonds  in 
the  structure  and  so  exerts  a  dominant  influence  on 
stability.  The  first-neighbor  constitution  also  sets  the 
local  atomic  packing  efficiency, which  controls  many 
physical  and  mechanical  properties,  including  the  mech¬ 
anisms  of  deformation.  The  first-neighbor  constitution  is 
given  by  partial  coordination  numbers  Zg,  which  are  the 
number  of  j  atoms  in  the  first  shell  of  a  reference  i  atom. 
In  binary  alloys  consisting  of  solute  (a)  and  solvent 
atoms  (Q),  there  are  four  separate  partial  coordination 
numbers:  Z^q,  Zq^,  and  Zqq.  Only  three  of  these  are 

independent,  because  Z^q  and  Zq^  are  related 

/aZaQ  =  /qZq,  [1] 

where  fi  is  the  atomic  fraction  of  atom  i.  Although 
recent  structural  descriptions  focus  on  solute-centered 
clusters,  a  complete  structural  description  requires  that 
equal  merit  be  given  to  both  solute-centered  and 
solvent-centered  clusters. 
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Partial  coordination  numbers  can  be  measured  using 
diffraction  techniques  (see  for  example.  Reference  3). 
Three  independent  measurements  are  needed  to  separate 
the  overlapping  contributions  to  the  total  scattering 
function  from  the  three  independent  atom  pairs  in 
binary  alloys,  and  six  independent  experiments  are 
needed  to  resolve  Zy  values  in  ternary  alloys.  An 
integrated  approach  combining  experimental  data, 
reverse  Monte  Carlo  calculations,  and  atomic  simula¬ 
tion  to  deconvolute  overlapping  peaks  has  been  devel¬ 
oped  recently  to  give  Zy  values  for  ternary  metallic 
glasses. The  Zy  values  for  more  complex  metallic 
glasses  have  been  estimated^^^  using  structural  concepts 
from  the  efficient  cluster  packing  (ECP)  model^^’^^  to 
rationalize  the  insignificance  of  some  weighting  factors. 
As  a  result  of  the  difficulties  associated  with  the 
determination  of  Zy,  relatively  little  data  are  available 
for  binary  glasses  and  almost  no  data  are  published  for 
partial  coordination  numbers  in  higher  order  glasses. 

Measured  partial  coordination  numbers  are  compiled 
in  a  recent  assessment  of  binary  metallic  glasses. These 
data  are  plotted  in  Eigure  1,  along  with  additional  data 
not  included  in  that  assessment. The  data  show 
significant  scatter,  and  a  critical  analysis  shows  two 
areas  of  concern.  Eirst,  several  datasets  do  not  satisfy 
the  required  relation  between  Z^q  and  Z^a  (Eq.  [1]).  In 
most  cases,  deviations  from  Eq.  [1]  are  probably  within 
experimental  error,  but  in  some  cases,  the  discrepancy 
can  be  much  larger  and  is  sometimes  as  large  as  3.  The 
second  issue  is  that  several  total  coordination  numbers 
around  a  and  Q,  Z^^  tot  and  Zq  tot?  deviate  significantly 
from  the  values  needed  for  efficient  atomic  packing  and 
consistency  with  measured  metallic  glass  densities.  In 
some  glasses,  the  total  coordination  number  is  as  low  as 
5  or  6,  once  again  suggesting  that  errors  in  some  Zy 
values  might  be  as  large  as  3.  The  experimental  data  are 
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Fig.  1 — Measured  partial  coordination  numbers  in  binary  metallic  glasses  as  a  function  of  the  atom  fraction  of  the  smaller  atom.  The  data  are 
grouped  by  the  relative  size  of  the  atoms.  Systems  with  a  nominal  radius  ratio,  R,  near  0.710  are  indicated  by  <9>  (the  maximum  number  of  lar¬ 
ger  atoms  that  can  be  first  neighbors  to  a  smaller  atom),  glasses  with  R  near  0.799  are  shown  as  <10>,  and  glasses  with  R  ^  0.902  are  shown  as 
<12>.  The  a  atom  is  taken  to  be  the  smaller  atom  and  is  the  larger  atom.  The  data  are  taken  from  Ref.  8  and  from  Refs.  5  and  9-11. 


given  in  Table  I,  and  the  values  with  high  assessed  errors 
are  indicated. 

In  addition  to  experimental  errors,  differences  in  Z^y 
could  come  from  systematic  structural  differences 
related  to  relative  atom  size  and  quench  rate.  The 
magnitude  of  chemical  interactions  between  the  atoms, 
Sfj,  might  also  influence  Zy  through  the  degree  of  short- 
range  atomic  order. Typically,  metallic  glasses  form 
in  systems  with  negative  heats  of  mixing,  suggesting  that 
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the  bond  energy  between  unlike  atoms  is  more  negative 
than  the  average  of  bond  energies  between  like-atom 
pairs,  eaQ<(^aa+  Sqq)/2.  Bonding  between  unlike  atoms 
is  increasingly  preferred  as  the  magnitude  of  this 
inequality  becomes  more  negative,  thus  influencing  Zy 
values.  The  largest  amount  of  partial  coordination  data 
is  available  for  Zr-Cu  and  Zr-Ni  systems.  Cu  and  Ni 
atoms  are  nearly  identical  in  size,  and  the  chemical 
interactions  in  these  two  binary  systems  are  expected  to 
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Table  I.  Experimental  Z^:  Data  are  Taken  from  Ref.  8  Unless  Otherwise  Noted 


Q 

a 

fa 

R 

Zaoi 

ZolQ. 

ZaoL 

Zq£1 

Z^a,tot 

Z^Q,tot 

Ref 

Ni 

B 

0.19 

0.698 

0 

9.3 

2.2 

10.8 

9.3 

13.0 

Ni 

B 

0.20 

0.698 

0 

5.8"^ 

1.5 

10.8 

5.8'^ 

12.3 

Ni 

B 

0.33 

0.698 

0.9 

9 

4.4 

9.4 

9.9 

13.8 

Ni 

B 

0.36 

0.698 

1.1 

8.7 

4.9 

9.2 

9.8 

14.1 

Fe 

B 

0.20 

0.704 

0 

8.6 

2.2 

12.4 

8.6 

14.6 

Y 

Cu 

0.33 

0.704 

2.9 

8.4 

4.1 

10.7 

11.3 

14.8 

Y 

Ni 

0.33 

0.704 

1.9 

9.3 

4.6 

9.5 

11.2 

14.1 

Y 

Cu 

0.833 

0.704 

10.4 

2.8 

14.0 

0 

13.2 

14.0 

Zr 

Be 

0.43 

0.709 

2 

lA 

5.4 

6.2 

9.1 

11.6 

Dy 

Ni 

0.30 

0.720 

1.3 

12 

3.1 

10.4 

8.5 

13.5 

Dy 

Ni 

0.30 

0.720 

10.8"^ 

4.6^ 

12.4"^ 

13.8'^ 

n.o"^ 

Dy 

Ni 

0.56 

0.720 

3.1 

5.2 

6.6 

6.9 

8.3 

13.5 

Ti 

Si 

0.16 

0.775 

0 

9.4 

1.8 

11.5 

9.4 

13.3 

Pd 

Si 

0.20 

0.775 

0 

6.6^ 

1.7 

10.6 

6.6^ 

12.3 

Y 

A1 

0.90 

0.788 

10.7 

1.6 

14.1 

1.1 

12.3 

15.2 

Ti 

Be 

0.38 

0.789 

2.9 

6.9* 

3.6* 

8 

9.8 

11.6 

Zr 

Ni 

0.33 

0.797 

1.3 

8.4 

4.2 

11.6 

9.7 

15.8 

Zr 

Cu 

0.35 

0.797 

4.6 

6.9 

3.7 

9.4 

11.5 

13.1 

10 

Zr 

Cu 

0.35 

0.797 

4.6 

lA 

4.0 

7.8 

12 

11.8 

Zr 

Cu 

0.35 

0.797 

2.8 

6.7* 

3.7* 

9.8 

9.5 

13.5 

Zr 

Ni 

0.35 

0.797 

3.3 

8.6* 

4.8* 

11 

11.9 

15.8 

Zr 

Ni 

0.36 

0.797 

2.3 

7.9* 

3.9* 

9.1 

10.2 

13 

Zr 

Ni 

0.36 

0.797 

3.3 

8.6 

4.8 

11 

11.9 

15.8 

Zr 

Cu 

0.40 

0.797 

2.0 

8.0 

5.4 

8.0 

10.1 

13.4 

5 

Zr 

Cu 

0.455 

0.797 

2.7 

lA 

6.2 

7.8 

10.1 

14.0 

5 

Zr 

Cu 

0.46 

0.797 

6.0 

5.5 

5.0^^ 

5^ 

11.5 

lO'^ 

Zr 

Cu 

0.46 

0.797 

3.2 

7.6 

6.5 

9.1 

10.8 

15.6 

9 

Zr 

Cu 

0.50 

0.797 

5.9 

6.0 

6.0 

8 

11.9 

14 

10 

Zr 

Cu 

0.50 

0.797 

5.8 

5.6 

5.6 

5 

11.4 

10.6 

Zr 

Cu 

0.50 

0.797 

4.9 

5.4 

5.4 

7.7 

10.3 

13.1 

Zr 

Cu 

0.50 

0.797 

3.0 

7.3 

7.3 

7.2 

10.3 

14.4 

5 

Zr 

Ni 

0.50 

0.797 

3.3 

6.7 

6.7 

7.8 

10 

14.5 

Zr 

Cu 

0.56 

0.797 

4.5 

6.5 

8.2 

6.1 

11.0 

14.3 

5 

Zr 

Cu 

0.57 

0.797 

5.4 

5.0* 

6.7* 

5.9 

10.4 

12.6 

Zr 

Cu 

0.614 

0.797 

6.1 

5.4 

8.6 

4.8 

11.5 

13.3 

5 

Zr 

Ni 

0.637 

0.797 

6.0 

5.0 

8.8 

5.8 

11 

14.6 

Zr 

Ni 

0.64 

0.797 

6.0 

5.0 

8.9 

5 

11 

13.9 

Zr 

Cu 

0.645 

0.797 

6.3 

5.5 

9.9 

4.9 

11.7 

14.8 

11 

Zr 

Cu 

0.645 

0.797 

7.2 

4.9 

9.0 

4.1 

12.1 

13.1 

5 

Zr 

Cu 

0.65 

0.797 

7.8 

4.6 

8.6 

6 

12.4 

14.6 

10 

Zr 

Cu 

0.65 

0.797 

6.9 

4.6 

8.6 

3.7 

11.5 

12.3 

Zr 

Cu 

0.65 

0.797 

6.7 

3.9* 

7.6* 

5.9 

10.6 

13.5 

Zr 

Ni 

0.67 

0.797 

6.4 

5.3 

10.8 

6.6 

11.7 

17.4 

Ni 

P 

0.20 

0.810 

0 

9.3 

2.3 

9.4 

9.3 

11.7 

Co 

P 

0.20 

0.816 

0 

8.9 

2.2 

10.1 

8.9 

12.3 

Fe 

P 

0.25 

0.816 

3.5 

8.1* 

2.6* 

10.4 

11.6 

13 

Pt 

Ge 

0.20 

0.820 

0 

5.6^ 

1.4 

10.3 

5.6'^ 

11.7 

Ta 

Ni 

0.50 

0.869 

4.9 

6 

6 

8.2 

10.9 

14.2 

Ta 

Ni 

0.55 

0.869 

6.4 

7.6* 

6.2* 

5.3 

14 

11.5 

Nb 

Ni 

0.40 

0.881 

3.8 

8.2 

5.5 

9 

12 

14.5 

Nb 

Ni 

0.50 

0.881 

5 

7.4 

7.4 

7.5 

12.4 

14.9 

Nb 

Ni 

0.56 

0.881 

5.5 

6.6 

8.4 

6.5 

12.1 

14.9 

Nb 

Ni 

0.62 

0.881 

6.1 

5.7 

9.3 

5.5 

11.8 

14.8 

Nb 

Ni 

0.63 

0.881 

6.6 

5.9 

10 

5.6 

12.5 

15.6 

Ti 

Cu 

0.50 

0.887 

4.5 

6 

6 

6.4 

10.5 

12.4 

Ti 

Ni 

0.40 

0.887 

2.3 

7.9 

5.3 

8.1 

10.2 

13.4 

Zr 

Pd 

0.30 

0.899 

3.6^^ 

5.4^ 

2.3 

11.5 

O.O'^ 

13.8 

*Data  do  not  satisfy  required  identity  in  Eq.  [1]. 

^Deviates  by  at  least  3  from  values  expected  for  efficient  topological  packing. 


be  roughly  similar  based  on  the  similar  electronegativ¬ 
ities,  compound  formation,  and  melting  temperatures  in 
these  two  systems.  Thus,  the  data  from  these  two 
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because  of  topology  and  chemical  effects.  Nevertheless, 
data  from  the  Zr-Ni  and  Zr-Cu  binary  systems  still 
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show  essentially  the  same  degree  of  scatter  (Figure  2)  as 
for  the  full  data  set  and  suggest  that  an  error  in  Z,y  of  ±1 
is  typical  and  may  be  as  large  as  ±1.5. 

Despite  significant  scatter,  clear  trends  are  observed. 
As  expected,  increasing  gives  higher  and  Zq^,  and 
decreases  Z^q  and  Zqq  systematically.  The  maximum 
values  of  and  Zqq  both  approach  approximately  12. 
Within  experimental  error,  Z^^  is  0  for  <  0.2  and  Zqq 
is  0  for  >  0.8  as  expected  in  systems  where  solute- 
solute  contact  is  avoided,  as  suggested  by  the  ECP 
model. The  maximum  values  for  Z^q  approach  ~9  to 
12  as  fa,  approaches  0,  and  the  maximum  in  Zq^ 


approaches  ~15  to  17  as  approaches  1.  Finally,  the 
total  coordination  around  a  ranges  from  approximately 
9  or  10  when  is  small  to  about  12  or  13  as 
approaches  1,  whereas  the  total  coordination  around  Q 
varies  from  approximately  12  or  13  when  is  small  to 
values  1 5  or  higher  as  approaches  1 . 

Although  these  trends  match  simple  expectations, 
currently  no  convenient  approach  is  available  to  esti¬ 
mate  Zij  in  metallic  glasses  that  accounts  for  relative 
atom  size  and  atom  fraction.  These  values  can  be 
determined  with  confidence  in  atomic  simulations,  but 
simulations  require  high-fidelity  atomic  potentials  that 
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Fig.  2 — Measured  partial  coordination  numbers  in  binary  Zr-Cu  and  Zr-Ni  metallic  glasses.  All  are  <10>  glasses  with  a  nominal  radius  ratio  of 
0.797.  The  data  are  taken  from  Ref.  8  and  from  Refs.  5  and  9-11. 
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are  difficult  to  produce  and  validate.  A  relative  lack  of 
high-quality  atomic  potentials  makes  it  difficult  to  use 
simulations  to  identify  trends  across  many  systems  with 
a  broad  range  of  topologies  and  chemistries.  A  recent 
effort  estimates  Z,y  in  binary  systems  by  including  effects 
of  chemical  short-range  ordering  into  cluster-packing 
models  via  a  quasi-chemical- type  approach.^  This 
gives  good  agreement  with  experimental  data  and  gives 
an  average  difference  from  measured  partial  coordina¬ 
tion  numbers  of  0.76,  which  is  a  slightly  better  than  an 
average  difference  of  0.82  obtained  from  a  simple 
topological  approach  that  does  not  include  chemical 
effects  (note  typographical  errors  in  the  text  and 
Table  III  of  Reference  12  regarding  ECP  predictions). 
However,  this  approach  requires  data  that  are  not 
widely  available.  The  ECP  model  gives  an  approach  for 
estimating  Zij  and  has  been  applied  successfully  to  many 
systems. However,  agreement  is  best  for  solute-lean 
glasses,  where  the  relative  atomic  size  and  solute 
concentrations  are  such  that  each  solute  can  have 
solvent  atoms  only  as  first  neighbors.  Simplifications 
in  the  approach  for  estimating  from  the  ECP  model 
gave  a  less  compelling  agreement  for  solute-rich  glasses, 
where  solute-solute  contact  is  unavoidable.  As  the  most 
stable  glasses  are  almost  always  solute  rich,^^^  extending 
a  predictive  capability  to  solute-rich  glasses  would  give 
an  important  new  capability. 

The  objective  of  the  current  work  is  to  develop  a 
topological  approach  for  estimating  partial  coordina¬ 
tion  numbers  in  binary  metallic  glasses  and  to  apply  this 
to  a  broad  range  of  binary  metallic  glasses.  An  explicit 
account  of  the  relative  atomic  size  and  atom  fraction  is 
included,  and  calculations  are  developed  over  the  full 
composition  range  so  that  solute-rich  glasses  are  cov¬ 
ered.  Predictions  are  compared  with  experimental  data. 
The  trends  in  Zy  as  a  function  of  relative  atomic  size  and 
composition  are  analyzed  to  explore  experimental 
scatter  and  to  look  for  new  insights  into  the  structure 
of  solute-rich  glasses. 


II.  DEFINITION  OF  STRUCTURAL  TERMS 


Structure-specific  equations  that  account  for  compo¬ 
sition  and  relative  atomic  size  are  developed  for  the  four 
Zij  values  in  binary  metallic  glasses.  Structural  terms 
that  include  topological  information  are  taken  from  the 
ECP  model.  ^  These  terms  are  introduced  briefly  and 
described  for  metallic  glasses  consisting  of  a  solute 
atoms  and  Q  solvent  atoms. 

The  ECP  model  allows  the  sites  in  a  metallic  glass 
structure  to  be  counted.  For  each  a  site  in  the  structure, 
there  is  1  ^  site,  2  y  sites,  and  Aq  Q  sites,  so  that  the  total 
number  of  structural  sites  per  a  site  is  (Aq  +  4) .  Sq 
ranges  from  8  to  20  and  depends  on  the  nominal  radius 
ratio  R  =  r^/ro^  of  solute  and  solvent  atoms. The  ECP 
model  also  allows  the  occupancy  of  sites  in  a  metallic 
glass  structure  to  be  specified.  Q  atoms  only  occupy  Q 
sites,  but  a  atoms  can  occupy  any  site  in  the  structure,  a 
and  Q  sites  are  always  occupied,  but  P  and  y  sites  can  be 
occupied  by  a  or  can  be  vacant.  The  a  atoms  first  occupy 


a  sites,  followed  by  P  sites,  then  y  sites,  and  finally  Q  sites 
as  the  solute  atom  fraction  increases. 

The  glass  is  solute  lean  when  the  number  of  a  atoms  is 
insufficient  to  fill  all  a,  P,  and  y  sites,  and  the  solute  atom 
fraction  is  given  as 

f.  =  S,/{Sa  +  S,)  [2a] 

where  Aa  is  the  total  number  of  a  atoms  in  the  struc¬ 
ture  normalized  by  the  number  of  a  sites.  In  solute- 
lean  glasses,  Aa  <  4.  Rearranging  terms  gives  Aa  as  a 
function  of  for  solute-lean  glasses 


In  solute-rich  glasses,  there  are  enough  a  atoms  to  fill 
all  a,  P  and  y  sites  and  some  of  the  Q  sites.  In  these 
structures. 


fa  =  S,/{Sa  +  ^) 

[3a] 

Sa=fa{Sa  +  ^) 

[3b] 

In  solute-lean  glasses,  Sq  is  constant.  However,  Sq 
changes  with  composition  in  solute-rich  glasses  as  Q 
sites  are  filled  progressively  by  a  atoms. If  a  is  smaller 
than  Q,  then  Sq  increases  with  increasing  and  Aq 
decreases  with  increasing/^^  when  a  is  larger  than  Q.  This 
is  accounted  for  through  an  effective  radius  of  atoms  on 
Q  sites,  tq 


where  S{ixq)  is  the  number  of  a  atoms  on  Q  sites 
normalized  by  the  number  of  a  sites.  In  this  weighted 
average,  (S{ixq)  / is  the  fraction  of  Q  sites  occupied 
by  a  and  {\  —  S(aQ) / Sq)  is  the  fraction  of  Q  sites 
occupied  by  Q.  S{aQ)  =  0  in  solute-lean  glasses  and 
A(aQ)  =  (Aa  —  4)  in  solute-rich  glasses.  The  influence  of 
atom  size  and  composition  is  obtained  by  using  an 
effective  radius  ratio  R  =  v^/vq  in  the  equation  for  Sq 
[see  Eq.  [7]  in  Reference  13].  Because  Sq  depends  on  R 
and  because  R  depends  on  Sq,  an  iterative  solution  for 
Sq  and  R  is  required. The  terms  R,  Sq,  S^,  and 
A(aQ)  are  evaluated  elsewhere  for  all  binary  metallic 
glasses. 

The  term  is  used  to  indicate  the  maximum 

number  of  Q  atoms  that  is  possible  in  the  first  shell  of 
another  Q  atom.  This  is  determined  by  setting  R  =  1  in 
the  equation  for  Sq,  which  gives  a  value  of  13.33.  This 
falls  between  the  maximum  number  of  unit  spheres  that 
can  simultaneously  contact  a  central  unit  sphere  (12) 
and  the  number  of  first-neighbors  in  a  body-centered 
cubic  structure  using  a  Voronoi  definition  of  coordina¬ 
tion  number  (14),  where  8  atoms  in  the  first  shell  contact 
the  central  atom,  and  the  remaining  6  atoms  are  at  a 
distance  approximately  1 5  pet  larger  than  the  minimum 
separation.  The  value  =  13.33  is  used  in  the 
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following  analyses  to  acknowledge  that  all  the  atoms  in 
the  first  coordination  shell  do  not  need  to  contact  the 
central  atom,  and  some  coordinating  atoms  may  be 
slightly  displaced  from  contact.  This  is  also  consistent 
with  the  approach  commonly  used  in  experiments, 
where  atoms  are  considered  as  first  neighbors  up  to 
the  first  minimum  in  the  radial  distribution  function, 
which  often  occurs  at  approximately  125  pet  of  the 
minimum  atomic  separation. 

The  terms  “a”  and  “solute”  have  thus  far  been  used 
interchangeably  and  are  taken  as  the  minority  species. 
However,  as  fa,  increases,  the  structure  will  eventually 
become  an  a-based  glass  with  Q  atoms  as  the  minority 
species.  This  transition  occurs  at  the  isostructural 
composition  fso^^^  When  the  a  and  Q  atoms  are  of 
equal  =  0.50.  However,  the  number  of  structural 

sites  in  metallic  glasses  depends  on  the  relative  sizes  of  a 
and  Q.  A  smaller  a  is  structurally  less  potent  because  it 
produces  fewer  Q  sites.  As  a  result,  f^^  is  greater  than 
0.50  when  R  <  1  and  is  less  than  0.50  when  R>  1.  In 
developing  equations  for  Zg,  the  terms  solute  and 
solvent  will  be  used  for  the  structural  minority  and 
majority  species,  respectively. 


III.  CALCULATION  OF  Zjj 

It  is  assumed  that  efficient  atomic  packing  occurs 
locally  around  both  solute  and  solvent  atoms,  as 
supported  by  phenomenological  correlations^^^  and  pre¬ 
vious  calculations. A  statistical  approach  is  used  to 
develop  the  Zg  equations  from  the  product  of  four 
terms:  (1)  the  number  of  reference  sites  occupied  by  an  i 
atom  normalized  by  the  number  of  a  sites  in  the 
structure,  (2)  the  probability  of  finding  an  i  atom  at  that 
site,  (3)  the  number  of  structural  sites  occupied  by  a  j 
atom  that  contact  the  i  site,  and  (4)  the  probability  of 
finding  a  j  atom  at  this  site.  This  product  gives  the  total 
number  of  i-j  contacts.  To  get  Z^,  the  number  of 
contacts  is  normalized  by  the  number  of  i  atoms  in  the 
system  considered.  A  more  detailed  description  of  the 
equations  is  given  in  the  Appendix. 

Partial  coordination  numbers  are  calculated  from 
Eqs.  [A4]  through  [A  10]  for  binary  systems  with  nom¬ 
inal  radius  ratios  of  7?  =  0.710,  0.799,  and  0.902.  These 
equations  are  applied  from  the  minimum  solute  con¬ 
centrations  needed  to  form  a  glass  (see  Eq.  [2a],  where 
Sa  =  1)  up  to  /iso  (0.60  for  R  =  0.710;  0.57  for 
R  =  0.799  and  0.53  for  R  =  0902).  When  /,  >/iso, 
the  calculations  are  performed  for  the  inverse  systems, 
where  a  becomes  the  solvent  and  Q  becomes  the  solute. 
Thus,  the  nominal  radius  ratio  in  the  inverse  system  R~^ 
is  given  by  1/7?  =  r^ijra,  and  the  solute  atom  fraction  in 
the  inverse  system  is  /“^  =  1  -/a.  Other  parameters  in 
Eqs.  [A4]  through  [A  10]  are  transformed  accordingly. 
To  avoid  confusion  and  to  enable  visualization  of  the 
trends  across  /so  while  maintaining  a  consistent  nota¬ 
tion,  the  experimental  Z^  values  in  Table  I  and  the 
calculated  Z^  values  in  the  following  figures  report  a  as 
the  smaller  atomic  species  over  the  full  range  of 

The  effective  radius  ratio  (ERR)  gives  a  second 
approach  to  estimate  Z^  based  on  the  assumption  that 
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a  global  alloy  composition  can  be  represented  by  the 
atomic  composition  of  an  /-centered  atomic  cluster,  with 
a  small  “stoichiometric  remainder”^^^  or  “glue  atoms”^^"^^ 
as  required.  Again,  efficient  local  atomic  packing  is 
assumed,  and  the  effective  radius  ratio  ^  =  r^/fQ 
described  previously  is  used  to  define  preferred  local 
cluster  coordination  at  integer  values  of  Z,  i,  and  / 
Equations  [A  12]  through  [A  17]  describe  the  ERR  cal¬ 
culation  to  generate  partial  and  total  coordination 
numbers  using  an  ideal  packing  efficiency  normalization 
and  linear  extrapolation  between  the  two  bounding, 
efficiently  packed  /-  or  /centered  clusters  as  a  function 
of /a. 


IV.  RESULTS  AND  DISCUSSION 

The  results  are  plotted  with  the  experimental  data  in 
Eigure  3.  Eor  clarity,  Zy  data  with  large  assessed  errors 
(Table  I)  are  omitted  from  this  comparison.  Predictions 
from  both  approaches  generally  match  the  experimental 
trends  and  reproduce  measured  data  within  the  assessed 
experimental  error  of  ±1.5.  The  calculated  values  satisfy 
Eq.  [1]  and  match  across  /so,  giving  additional  confi¬ 
dence  in  the  estimates.  Chemical  effects  are  expected  to 
influence  the  local  atomic  structure,  and  detailed 
approaches  such  as  that  of  Lass  et  are  expected 
to  be  useful.  However,  the  current  results  show  that 
topology  alone  gives  a  good  representation  of  the  data, 
suggesting  that  the  magnitude  of  chemical  effects  on  Zy 
values  may  be  within  experimental  scatter. 

The  development  of  equations  for  Zy  from  the  ECP 
model  uses  a  statistical  view  of  topology,  and  the  Zy 
values  calculated  here  are  global  averages.  In  the  ECP 
approach,  the  same  type  of  atom  can  have  different  Zy 
values  when  occupying  different  sites.  Eor  example,  in 
the  same  structure,  the  number  of  a  atoms  around  an  a 
atom  that  occupies  an  a  site  can  be  different  than  the 
number  of  a  atoms  around  an  a  atom  that  occupies  an  Q 
site.  The  differences  are  not  large.  Eor  example,  in 
solute-rich  glasses  with  a  nominal  R  of  0.799  and 
^(o^q)  ^  2,  Zaoi  is  2  for  a  atoms  on  a  sites  and  is  1.9  for  a 
atoms  on  Q  sites.  At  a  higher  so  that  A(aQ)  =  7,  Z^a, 
is  7  for  a  on  a  sites  and  is  5.9  for  a  on  Q  sites.  Similar 
differences  are  found  for  Z^q  and  Zq^  in  solute-rich 
glasses.  No  other  site-specific  differences  in  Zy  are 
indicated  by  Eqs.  [A4]  through  [A  10]  in  the  Appendix. 

The  ERR  approach  gives  cluster  compositions  that 
correspond  closely  to  the  bulk  alloy  composition,  where 
efficient  local  atomic  packing  is  obtained  with  an 
integral  number  of  constituent  atoms  in  the  first 
coordination  shell.  A  strong  correlation  is  found 
between  these  compositions  and  glass-forming  ability, 
suggesting  that  these  efficiently  packed  /-  or  j-  centered 
clusters  may  be  statistical  representatives  of  majority 
structural  elements  of  the  glass. Although  the  ERR 
approach  implicitly  assumes  that  the  first-neighbor 
environments  around  a  given  species  are  equivalent,  it 
remains  to  be  shown  that,  for  example,  the  coordination 
shell  of  an  efficiently  packed,  a-centered  cluster  is 
identical  to  the  first-neighbor  environment  of  an  a  atom 
in  the  first  shell  of  such  a  cluster.  Even  if  this  is  found  to 

METALLURGICAL  AND  MATERIALS  TRANSACTIONS  A 


6 

Distribution  Statement  A.  Approved  for  public  release;  distribution  unlimited. 


20 

15 

5 

0 

0.0  0.2  0.4  0.6  0.8  1.0 

Atom  fraction  a,  fa 


(b) 


(e) 


0.0  0.2  0.4  0.6  0.8  1.0 


Atom  fraction  a,  fa 

(c) 


Atom  fraction  a,  fa 

(0 


Fig.  3 — Calculated  partial  coordination  numbers  in  binary  metallic  glasses  compared  with  experimental  data.  The  experimental  data  are  taken 
from  Ref.  8  and  from  Refs.  5  and  9-11. 


be  true,  the  condition  where  first-neighbor  environments 
are  identical  for  a  given  species  throughout  the  structure 
applies  only  to  discrete  compositions.  For  most  com¬ 
positions,  the  structure  consists  of  efficiently  packed 
clusters  with  a  stoichiometric  remainder  of  atoms  or  a 
mixture  of  more  than  one  type  of  efficiently  packed 
clusters.  Currently,  there  is  no  way  to  estimate  the  first- 
neighbor  environment  of  the  stoichiometric  remainder, 
and  there  is  no  reason  to  expect  that  the  first  shell  of 
these  atoms  will  be  identical  to  the  other  atoms.  Thus,  a 
given  glass  structure  may  have  a  range  of  first-neighbor 
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environments  with  local  Z,y  values  that  vary  slightly 
from  the  global  average  Zfj  values. 

From  the  previous  discussion,  Zg  values  from  both 
the  ECP  and  ERR  approaches  suggest  that  small, 
local  variations  in  first-neighbor  environments  are 
expected.  This  is  consistent  with  the  concept  of 
quasi-equivalent  clusters,^^^^  where  distributions  in 
first-neighbor  environments  are  shown  in  atomistic 
simulations.  This  gives  topological  fluctuations  in  both 
local  constitution  and  packing  that  are  expected  to  be 
important. 


VOLUME  43A,  AUGUST  2012—2655 


7 

Distribution  Statement  A.  Approved  for  public  release;  distribution  unlimited. 


Predictions  for  Zq^  from  the  ECP  model  diverge  from 
experimental  data  when  R  <  I  and  >  /iso-  This  repre¬ 
sents  structures  where  the  larger  Q  atoms  are  surrounded 
increasingly  by  the  smaller  a  atoms.  These  structures  can 
be  viewed  equivalently  as  solute-lean  glasses  with  R  >  1, 
where  larger  solute  atoms  are  surrounded  by  smaller 
solvent  atoms.  Like  glasses  with  R  <  1,  these  solutes  are 
expected  to  occupy  a,  and  y  sites  progressively  as  the 
atom  fraction  of  the  larger  solute  atoms  increases.  The 
smaller  experimental  Zq^  values  suggest  that  contact 
between  larger  atoms  occurs  at  lower  concentrations 
than  predicted  by  the  ECP  model,  so  that  the  larger 
atoms  might  not  occupy  ^  and  y  sites.  To  explore  this 
possibility,  the  sizes  of  P  and  y  sites  are  estimated  in  the 
Appendix  as  a  function  of  the  solute  size  and  are  plotted 
in  Figure  4(a).  Although  the  sizes  of  P  and  y  sites  show 
significant  variability  because  of  the  nonspherical  nature 
of  the  clusters,  in  general  the  P  sites  can  accommodate  a 
atoms  up  to  7?  1.12.  When  solute  atoms  are  larger  than 

this,  the  a  atoms  can  no  longer  be  placed  in  P  sites 
without  pushing  apart  the  constituent  clusters.  This  is 
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Fig.  4 — (a)  Estimated  size  (see  Appendix)  of  atoms  that  can  fit  into 
P  and  y  sites  without  forcing  clusters  apart.  The  lower  bounds  for 
the  P  site  and  y  site  bands  are  for  face  contact,  the  upper  bounds  are 
for  edge  contact,  and  the  centerline  is  a  simple  average  of  these  two 
bounds.  The  dotted  line  represents  the  condition  where  the  site  size 
is  exactly  equal  to  the  size  of  the  a  atom.  and  the  solute  site  sizes 
are  both  relative  to  the  radius  of  Q  atoms,  tq.  (b)  Illustration  show¬ 
ing  the  size  of  P  sites  relative  to  the  size  of  the  solutes.  When  <  1, 
the  solute  can  fit  into  P  sites,  but  when  >  1,  the  solute  is  too  large 
to  fit  in  p  sites  without  pushing  apart  the  bounding  clusters. 


illustrated  in  Figure  4(b),  where  a  solute  with  R  <  I  fits 
in  a  ^  site  but  a  solute  with  R>  I  overlaps  atoms  in  the 
bounding  clusters.  In  contrast,  the  y  sites  are  too  small  to 
hold  a  atoms  without  increasing  the  cluster-cluster 
separation  over  the  full  range  of  a  atom  sizes.  Figure  4(a) 
suggests  that  filling  y  sites  with  a  solutes  will  produce 
increasingly  larger  strains  between  structure-forming 
clusters  with  increasing  R. 

Experimental  evidence  gives  clues  to  the  site-filling 
sequence  in  metallic  glasses  and  supports  this  suggestion 
that  a  occupancy  of  y  sites  increases  cluster-cluster 
separation.  Diffraction  results  show  that  a-a  contact 
does  not  occur  in  any  glass  when  R  <  1  and  <  0.30 
(Table  I),  demonstrating  a  remarkable  tendency  for 
solute  atoms  to  avoid  contact.  In  the  ECP  model,  this 
avoidance  is  accomplished  by  filling  a,  P,  and  y  sites 
between  solute-centered  clusters  with  solvent  atoms  only 
in  the  first  shells.  Because  Q  sites  surround  a,  P,  and  y  sites, 
solute-solute  contact  can  only  occur  when  oc  atoms 
occupy  Q  sites.  Solute-solute  contact  becomes  unavoid¬ 
able  when  R  <  I  and >  0.30,  after  the  a,  P,  and  y  sites 
are  fully  occupied  and  additional  solute  atoms  must 
occupy  Q  sites.  Consistent  with  this,  the  data  in  Table  I 
show  that  solute-solute  contact  occurs  in  every  glass 
where  /a  >  0.30.  Additional  support  for  the  filling  and 
relative  sizes  of  P  and  y  sites  is  given  by  density 
measurements  in  melt-quenched  Fe-B  metallic  glass 
ribbons.  The  density  is  relatively  insensitive  to  concen¬ 
tration  from  boron  atom  fractions  of  0.13  to  approxi¬ 
mately  0.19,  and  the  density  drops  linearly  from  boron 
atom  fractions  of  approximatel^y  0.19  to  0.27.^^^^  In  Fe-B, 
the  P  sites  are  fully  occupied  {S^  =  2)  between  =  0.18 
and  =  0.19  (Table  A1  in  Reference  8),  and  the 
decreasing  density  at  higher  boron  atom  fractions  is 
consistent  with  the  idea  that  the  structure-forming  clus¬ 
ters  are  pushed  apart  as  soon  as  boron  begins  to  fill  y  sites. 

The  most  stable  Fe-B  glasses  have  a  boron  atom 
fraction  in  the  range  of  0.18  to  0.20,  so  that  the  internal 
strains  and  free  volume  associated  with  filling  of  y  sites 
may  contribute  to  the  decreased  stability.  Nevertheless, 
nearly  three  quarters  of  the  y  sites  are  filled  by  boron 
before  Fe-B  glasses  can  no  longer  be  made.  These 
observations  support  the  concept  that  a,  P,  y,  and  Q  sites 
are  filled  progressively  by  solutes  with  increasing  in 
structures  with  R  <  1,  and  that  solute  occupancy  of  P 
sites  does  not  distort  the  structure  while  solute  occu¬ 
pancy  of  y  sites  introduces  a  structural  expansion. 

Only  two  Zij  data  sets  are  available  for  solute-lean 
glasses  with  R>  I  (see  Table  I  for  glasses  with  R  <  I 
and  /a  >  0.79),  and  the  values  are  ambiguous  regarding 
avoidance  of  solute-solute  contact.  Y-Y  contact  is 
absent  in  Yo.i67Cuo.833  at  a  composition  that  suggests 
that  Y  fills  a,  P,  and  y  sites,  whereas  Zyy  =  IT  in 
Yo.ioAlo.9o  at  a  composition  where  the  Y  fills  only  a  and 
P  sites.  Errors  in  Zyy  from  Y0.10AI0.90  are  estimated  to 
be  ±50  pet  as  a  result  of  counting  statistics, and 
additional  errors  are  expected  from  the  assumed  atomic 
radii  used  to  fit  the  data.^T  No  data  are  available  for 
glasses  with  R  >  1  between  the  solute  atom  fractions  of 
0.167  < /a  <  0.33  (see  Table  I  for  glasses  with  R  <  I  and 
0.67  < /a  <  0.833).  However,  the  site-filling  sequence  can 
be  inferred  from  Figure  5,  which  shows  that  most 
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Fig.  5 — Topologies  (relative  sizes  and  concentrations)  for  all  binary 
metallic  glasses  showing  a  preference  for  metallic  glasses  with  R  <  1. 
Most  binary  glasses  with  R>  1 .26  have  solute  atom  fractions  sug¬ 
gesting  that  a  and  P  sites,  but  not  y  sites,  are  occupied  by  solutes 
(Soc  <  2).  Redrawn  from  Figure  5  in  Ref.  8. 

known  binary  metallic  glasses  have  R  <  1.  These  glasses 
occur  over  the  full  range  of  solute  atom  fractions,  and 
the  discussion  supports  the  view  that  a,  1^,  y,  and  Q  sites 
are  filled  progressively  as  increases.  Binary  glasses 
with  \<R<\.26  are  also  often  solute  rich  (5a  >4), 
suggesting  that  solutes  occupy  a,  ^5,  y,  and  Q  sites  in 
these  glasses  as  well.  In  contrast,  almost  all  binary 
glasses  with  R  >1.26  are  solute  lean,  with  solute  atom 
fractions  that  allow  only  a  and  P  sites  to  be  filled  by 
solutes  (5a  <  2).  This  suggests  that  only  a  and  P  sites  are 
filled  in  binary  glasses  with  R  >  1.26,  and  that  expan¬ 
sion  caused  by  placing  larger  solutes  on  y  sites  destabi¬ 
lizes  the  amorphous  structure.  This  gives  a  structural 
asymmetry  that  can  explain  the  observed  preference  for 
binary  glasses  with  R  <  I  and  represents  an  important 
modification  to  the  ECP  model. 

This  change  in  site  filling  can  be  modeled  by  a  simple 
modification  in  the  equations  used  to  predict  Z,y.  Because 
the  divergence  applies  to  values  of  Zoa  in  glasses  where 
R  <  I  and  /a  >  /iso,  the  equations  for  ZaQ  are  modified 
for  glasses  where/,  < /so  and  R>1.  Equation  [A8]  is 
unchanged  but  is  valid  only  to  <2  to  reflect  the 
change  in  the  boundary  between  solute-lean  and  solute- 
rich  structures.  This  same  change  is  applied  in  Eq.  [A9] 
by  replacing  all  of  the  4s  with  2s.  The  predictions  from 
this  simple  modification  show  better  agreement  with 
experiment  (compare  Figure  6  with  Figure  3(d)). 

From  the  ERR  approach,  the  structure  can  be  defined 
also  by  efficiently  packed  clusters  with  mixed  first-shell 
occupancy  and  only  a  relatively  small  number  of 
additional  atoms  as  a  stoichiometric  remainder  between 
these  clusters.  Currently,  there  is  no  analysis  for 
determining  whether  the  ERR  and  ECP  structures  are 
essentially  identical  or  for  distinguishing  between  these 
possibilities.  Although  based  on  ECP  concepts,  the  ERR 
approach  differs  in  two  major  ways:  (1)  It  allows  solute- 
solute  contact  by  assigning  the  first  addition  of  a  solute 
atom  beyond  the  minimum  concentration  needed  to 
form  solute-centered  clusters  to  the  first  shell  of  the 
structure-forming  clusters,  and  (2)  it  does  not  account 
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Fig.  6 — Calculated  Zq^^,  modified  to  account  for  the  inability  to  fill  y 
sites  in  glasses  with  R  <  1  and  A  >  /iso,  compared  with  experimental 
data.  The  experimental  data  are  taken  from  Ref.  8  and  from  Refs.  5 
and  9-11. 


specifically  for  intercluster  sites  such  as  P  and  y  sites.  A 
shortcoming  of  this  technique  is  that  at  lower  solute 
contents  where,  by  thermodynamics,  solute-solute  con¬ 
tact  is  avoided,  the  ERR  method  overpredicts  the  first 
neighbor  solute-solute  coordination  contribution.  Nei¬ 
ther  the  ERR  nor  the  ECP  models  specify  the  arrange¬ 
ment  of  atoms  within  the  coordination  shell. 

The  current  equations  predict  the  partial  coordination 
numbers  to  within  experimental  scatter,  but  they  are 
rather  cumbersome  and  require  additional  calculations 
to  convert  /,  to  structure-specific  ECP  parameters  or 
require  the  consideration  of  many  permutations  of 
possible  atomic  clusters.  As  an  aid,  simple  phenomeno¬ 
logical  equations  to  estimate  the  partial  coordination 
numbers  of  binary  metallic  glasses  as  a  function  of 
composition  and  relative  atom  size  are  offered  in  this 
study.  As  a  function  of /,  only,  the  equations  are 

Zaa  =  14.8/a  “  2.8  [5a] 

ZaQ  =  11.3 -9.8/a  [5b] 

ZQa  =  16.8/a -1.7  [5c] 

Zqq=  14.2-14.1/a  [5d] 

Both  /a  and  R  are  accounted  for  independently  in  the 
phenomenological  equations 

Zaa  =  14.57/a  +  2.79R  -  4.91  [6a] 

ZaQ  =  -9.91/a  +  1.37R  +  10.26  [6b] 

ZQa  =  16.90/a  -  1.47R  -  0.53  [6c] 

Zqq  =  -14.28/a  +  2.32R  +  12.43  [6d] 
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These  simple  phenomenological  equations  predict 
within  experimental  scatter,  and  the  squared  linear 
regression  correlation  coefficients  range  from  0.85  to 
0.94  for  all  equations.  The  equations  for  and  Zqq  are 
undefined  for  values  less  than  zero,  and  the  equations  for 
Z^Q  and  Zq^  are  undefined  for  values  less  than  one. 


V.  CONCLUSIONS 

The  equations  for  partial  coordination  numbers  in 
binary  metallic  glasses  are  derived  from  statistical 
topology  and  the  requirement  of  efficient  local  atomic 
packing.  Two  related  topological  models  are  used,  the 
ECP  model  and  the  ERR  model.  The  equations  from 
both  ECP  and  ERR  topologies  give  a  good  ability  to 
estimate  systematic  trends  in  Zij  as  a  function  of  relative 
atomic  size  and  composition.  The  estimates  are  typically 
accurate  to  within  the  assessed  experimental  scatter  of 
±1.5,  validating  the  underlying  structural  concepts. 
Although  expected  to  be  important,  the  magnitude  of 
a  chemical  influence  on  local  atomic  structure  is  less 
than  the  current  experimental  scatter.  Deviations  from 
experimental  data  suggest  improvements  in  the  ECP  and 
ERR  structural  topologies.  The  ECP  model  overesti¬ 
mates  Zq^  when  the  atom  fraction  of  a  approaches  1, 
leading  to  the  observation  that  P  sites  can  accommodate 
a  solutes  up  to  R  <  1.12,  but  an  expansion  of  bounding 
clusters  is  needed  to  place  a  solutes  in  P  sites  when 
R>  1.12.  y  sites  are  always  smaller  than  oc  solutes,  so 
that  it  is  increasingly  difficult  to  place  a  solutes  in  y  sites 
as  solute  size  increases.  Observed  glass  topologies  and 
estimates  of  the  sizes  of  P  and  y  structural  sites  in  the 
ECP  model  suggest  that  y  sites  are  not  occupied  when 
R  >  1.26.  Modified  equations  give  an  improved  match 
with  experimental  data  when  this  topological  anomaly  is 
considered.  This  important  modification  gives  an  asym¬ 
metry  to  the  site-filling  rules  of  the  ECP  model  and 
explains  a  preference  for  glasses  with  R  <  1.  The  current 
application  of  the  ERR  model  does  not  account  for 
overlapping  structure-forming  clusters,  and  hence,  it 
does  not  accommodate  the  experimentally  observed 
tendency  to  avoid  solute-solute  contact  for  solute-lean 
metallic  glasses.  Simplified  empirical  equations  for  Zij 
are  given  for  convenience. 
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APPENDIX 

An  accounting  of  all  possible  contacts  in  a  glass 
system  can  be  summarized  by  the  equation 


S1P1S2P2  ±  S2P3S4P4  ±  SspsSePe 


[Al] 


where  and  p^  are  the  number  of  reference  sites  and 
their  probability  of  occupancy,  whereas  m  is  the 
number  of  reference  atoms  of  interest.  For  consistency 
across  Zg,  all  refer  to  the  same  site  in  all  cases;  Si 
describes  the  occupancy  of  a  sites,  whereas  S2  describes 
the  species  occupying  the  Q  sites  coordinated  around 
that  a  site.  A3  describes  the  occupancy  of  the  Q  sites, 
whereas  S4  describes  the  occupancy  of  a  sites  coordi¬ 
nated  around  that  Q.  Likewise,  A5  describes  the 
occupancy  of  the  Q  sites,  whereas  describes  the 
occupancy  of  Q  sites  around  an  Q.  Using  this 
approach,  two  relationships  are  developed  for  each 
Zif.  one  for  solute-lean  glasses  (Aa  <  4)  and  one  for 
solute-rich  glasses  (Aa>4).  Table  AI  lists  the  values  of 
Sn  and  Pn  for  each  case  of  interest. 

As  an  example,  consider  Z^^  for  the  solute-rich  case 
where  a-a  contact  can  occur  three  different  ways:  from 
an  a  atom  on  a  solute  site  to  an  a  atom  on  an  Q  site, 
from  an  a  atom  on  an  Q  site  to  an  a  atom  on  a  solute 
site,  and  from  one  a  atom  on  an  Q  site  to  another  a 
atom  on  a  different  Q  site.  Because  solute  sites  are 
always  separated  by  Q  sites,  there  is  no  a-a  contact 
between  solute  sites.  Starting  with  the  first  of  these 
terms,  we  define  values  for  each  of  the  four  terms  (1) 
through  (4)  described  previously  that  correspond  to 
S\,P\,S2,  and  P2,  respectively.  On  average,  four  solute 
sites  contact  each  Q  site  (Ai),  and  the  probability  of 
these  sites  being  occupied  by  an  a  atom  is  1  {pi).  The 
number  of  Q  sites  around  each  solute  site  (A2)  is  Aq, 
and  the  probability  of  finding  an  a  atom  on  one  of 
these  sites  {p2)  is  S{(Xa)/Sa.  The  product  of  these  four 
terms  is 


S1P1S2P2  =  (4)(1)(Aq)  ^  ^  4A(aQ)  [A2] 


The  same  result  is  obtained  for  the  number  of  a-a 
contacts  that  originate  from  an  a  atom  on  an  Q  site  to 
an  a  atom  on  a  solute  site  (A3J93A4/74).  To  determine  the 
third  product,  there  are  Aq  Q  sites  per  a  site  (A5),  and 
the  probability  of  an  a  atom  on  one  of  these  sites  (ps) 
is  S{(Xq)/Sq.  To  estimate  the  number  of  Q  sites  that 
contact  each  Q  site,  we  note  that  on  average  four  solute 
sites  contact  each  Q  site,  and  these  sites  are  all 
occupied  by  a  atoms  in  solute-rich  glasses.  The  total 
number  of  a  atoms  that  can  pack  around  an  Q  atom  is 
given  by  Aq  evaluated  at  fQ/r^  =  1/R,  which  is  denoted 
as  Aq^.  The  fraction  of  space  in  the  first  shell  of  each 

Q  atom  occupied  by  these  four  solutes  is  4/Aq^,  and 
the  fractional  space  remaining  for  Q  atoms  is 
1— 4/Aq^.  The  total  number  of  Q  atoms  that  can 
coordinate  with  a  given  Q  atom  (A6)  is  thus 

(1  “  4/^q^)-  The  probability  of  an  a  atom  on 

one  of  these  sites  (p^)  is  S{(Xq)/Sq.  The  product  of 
these  four  terms  is 
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Table  AI.  Values  of  Reference  Sites  (SJ  and  Occupation  Probabilities  (pj 


S, 

5, 

Pi 

52 

P2 

53 

P3 

54 

P4 

55 

Ps 

a  on 

a  on 

a  on 

a  on 

a  on 

a-site 

Q-site 

Q-site 

a-site 

Q-site 

4  4 

^  4 

Sn 

0 

5n 

0 

4  4 

^  4 

Sn 

0 

>4 

4 

1 

Sn 

5'(an) 

Sn  ^ 

Sg 

4 

1 

Sn 

5'(aa) 

4 

Pe 


<4 

>4 


a  on 
a-site 


<4 

4  1 

5n 

>4 

4 

1 

Sn 

ZjJa 

<4 

a  on 

a-site 

4  f 

Sn 

>4 

4 

1 

Sn 

Q  on 
Q-site 
1 


_  >^(an) 
Sq 


Q  on 
Q-site 
1 


_  SjoCg) 
Sq 


Sn 

Sa 

Sa 

Sa 

Sn 


Q  on 
Q-site 
1 


1  _ 


a  on 
Q-site 
0 

Sjocg) 

Sg 

Q  on 
Q-site 
1 


0^=1 


0^=1 


0^=1 


0^=1 


0^=1 


0^=1 


&R=l 


Sn  1  -  ^ 


^R=l 


a  on  Q-site 
around  Q-site 
0 


=1  1  _ 

1  M/R 


1 

-I  /R 


Sjocg) 

Sg 


Q  on  Q-site 
around  Q-site 

1  ' 

^  ^1/7? 


1-4^ 


1  _ 


Q  on  Q-site 
around  Q-site 

1-'  ■  ' 


1-^ 


1 

Sg 


(X  on  Q-site 
around  Q-site 
0 


1 

^  -1  /R 


Sjocg) 

Sg 


#  a-sites 

4 

S{xq)  +  4 

#  Q-sites 
Aq 

Aq  -  A(aQ) 

#  a-sites 
Aa 

A(aQ) +  4 

#  Q-sites 
Aq 

Sq  -  S{(Xq) 


*When  <  1  and  >fso^  all  values  of  4  are  replaced  by  values  of  2. 


SsPsSePe  =  (^q)  ^ 


[A3] 


¥) 


^Soc<4 

^aQ 


[A8] 


Finally,  we  add  the  three  previous  products  and 
normalize  by  the  number  of  a  atoms  per  a  site  in  the 
structure  to  give  the  solute-rich  a-a  partial  coordination 
number 


'^Soi  >  4 
^aa 


45(an)  +  45(o:n)  +  (5(aa))"  (^) 
S{Gia)  +  4 


[A4] 


Similar  analyses  result  in  the  following  complete  set  of 
expressions  for  other  partial  coordination  numbers  for 
solute-rich  and  solute-lean  glasses  that  also  can  be 
shown  to  satisfy  Eq.  [1] 

=  0  [A5] 

=  |A61 


ySa  >4  _ 

^aQ  ~ 


(^l_m){4Sa  +  S{an) 


Van) +4 


[A9] 


=  V  [AlO] 

Equations  [A8]  and  [A9]  are  linked  to  Eqs.  [AlO]  and 
[All]  through  the  relation  in  Eq.  [1].  It  can  be  shown 
that  fa/ fa  =  S^/Sq  for  solute-lean  glasses,  so  that 
Eq.  [1]  is  satisfied  for  Eqs.  [A8]  and  [AlO]  by  inspection. 
A  more  detailed  comparison  shows  that  Eqs.  [A9]  and 
[All]  also  satisfy  Eq.  [1]. 

The  ERR  approach  uses  specific  a-  and  Q-centered 
clusters  that  define  the  base  alloy  composition,  where 
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mixed  occupancy  in  the  coordinating  shell  is  always 
acceptable.  This  methodology  essentially  deals  with 
integer  total  and  partial  coordination  numbers  for  each 
species  around  a  specified  central  atom  and  can  be 
extended  easily  into  multicomponent  alloy  systems. 
Similar  to  R  =  r^/r^  where  the  radius  of  a  is  smaller 
than  that  of  Q,  the  effective  radius  ratio  for  an 
a-centered  cluster  {R)  is  defined  by: 


^  _  R{^aa  + 

Z^^R  +  ZaQ 


[A12] 


Equation  [A  12]  must  also  satisfy  local  coordination 
summation  where  =  Z^^  tot-  Here  Z^^  tot  is  an 

integer  corresponding  to  the  total  coordination  number 
around  the  a  atom  in  an  a-centered  cluster.  Further¬ 
more,  the  partial  coordination  numbers  Z^a  and  Z^q  can 
be  described  in  terms  of  atom  fraction  of  a  (/«)  in  an 
a-centered  cluster  by  Zaa  =  1/a  (Za, tot  +  l)~l|  and 

ZaQ  —  Zq(  tot  l/a  (■Z^a,tot  1 )  1 1  ?  hcnCC 

^  ^ _ ^(^Mot) _ 

“  /a(Za,tot  +  1)-1||^+  [Za,tot  -  /a  (Za,tot  +  l)-l|] 

[A13] 


where  by  mathematical  limits,  if  R  <\,  then  0  < 
tot  <  13  and  14  <  za,tot  ^20.  It  is  assumcd  that  the  total 
and  partial  coordination  numbers  of  each  species  Z^^  tot, 
-^Q,tot,  -^aa,  -^aQ,  and  Zqq  adjust  linearly  as  a 

function  of  fa,  between  integer  values  of  Z^^  tot  and  Zq  tot- 
So  it  is,  through  Z^^  tot  and  Zq  tot,  thus  possible  to 
estimate  total  and  partial  coordination  numbers  as  a 
function  of for  a  given  radius  ratio  R,  which  is  shown 
in  Figure  3  for  R  values  of  0.902,  0.799,  and  0.7099. 

To  estimate  the  sizes  of  P  sites  (interstices  at  the  center 
of  an  octahedron  of  structure-forming  clusters)  and  y 
sites  (interstices  at  the  center  of  a  tetrahedron  of 
structure-forming  clusters),  we  start  by  approximating 
the  cluster  radius  as  one  half  the  cluster-cluster  separa¬ 
tion.  The  minimum  cluster-cluster  separation  represents 
the  most  efficient  packing  and  occurs  when  the  faces  of 
adjacent  clusters  are  in  contact.  Face  contact  may  not  be 
possible  for  all  cluster  contacts,  and  so  an  intermediate 
separation  for  edge  contact  of  adjacent  clusters  is  also 
considered.  The  equations  for  these  cluster-cluster 
separations  are  taken  from  Eqs.  [A3 (a)]  and  [A3(d)]  in 
Reference  7  and  after  normalization  by  tq  are 

^dSter  =  (1 


Similarly,  for  Q-centered  clusters  where  the  radius  of  a 
is  smaller  than  that  of  Q 


^  _  (Zpa  +  Zqq) 

^  ZaaR  +  Zqq 


[A14] 


where  Zq^  +  Zqq  =  ZQ^tot  and  ZQ^tot  is  an  integer  cor¬ 
responding  to  the  total  coordination  number  around 
an  Q  atom  in  an  Q-centered  cluster.  The  partial  coor¬ 
dination  numbers  Zq^  and  Zqq  can  be  described  in 
terms  of  in  an  Q-centered  cluster  by 

Zaa  —  1/a  (Za, tot  T  l)  1|  and  Zqq  —  ZQ  ^ot  /a(ZQ^tot“t“ 
1);  hence 

^  ^  _ [Zg^iOi _ 

""  [/a(Za,tot  +  1)]^  +  [ZQ,tot  -/a(ZQ,tot  +  l)] 

[A15] 


<?ter  =  [A19] 

The  distance  from  the  center  of  a  site  to  the  center 
of  an  a  site  of  one  of  the  bounding  clusters  (Figure  4(b)) 
is  simply  \/2i?ciuster  -  This  same  distance  is  also  given  by 
the  distance  between  an  a  and  P  atom  separated  by  a 
trigonal  ring  of  Q  atoms  (a  trigonal  bipyramid)  for  face 
contact,  and  for  the  distance  between  an  a  and  P  atom 
separated  by  a  pair  of  contacting  Q  atoms  for  edge 
contact.  These  latter  two  relations  are  taken  from 
Eqs.  [A  1(a)]  and  [A  1(d)]  in  Reference  7. 

<T  =  [A20] 


Rarely  will  R  be  equal  to  the  ideal  radius  ratio  for 
perfect  packing  for  a  given  coordination  number,  so  to 
normalize  coordination  numbers  to  perfect  packing,  R  is 
equated  to  the  R  value  for  a  specific  a  total  coordina¬ 
tion  number,  whereby  a  specific  atomic  fraction  of  a  can 
be  calculated  for  each  total  coordination  number 
(integer)  for  a  given  radius  ratio  R.  That  is. 


and 


fa 


/?(^+l)  -(Z,,tot+l) 
(^-  l)(Z,,,o,  +  l) 


(R  -  l)(Za,tot) 


[A16] 


[A17] 


—  fJZZf-l  +  \/(l  +  Rfi)  -1  [A21] 


Setting  these  equations  equal  to  i/2-Rciuster  and  solving 
for  Rp  gives 


oface  _ 

Rp  - 

[A22] 


^edge  ^ 


(V2-i)\/(1+/?„)2-1  +  1 


+  l  -  I 

[A23] 
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A  similar  approach  is  used  to  estimate  the  size  of 
atoms  that  can  occupy  y  sites.  In  this  case,  where  a  y  site 
is  surrounded  by  a  tetrahedron  of  solute-centered 
clusters,  the  distance  from  the  center  of  a  7  site  to  the 
center  of  an  a  site  of  one  of  the  bounding  clusters  is 

Y^^ciuster.  The  development  is  otherwise  the  same  as  for 
^  sites,  and  the  final  equations  are 


[A25] 


Equations  [All]  and  [A23]  are  plotted  in  Figure  4(a) 
as  the  upper  and  lower  bounds  for  the  size  of  atoms  that 
can  be  placed  at  P  sites  without  pushing  apart  the 
bounding  clusters,  and  Eqs.  [A24]  and  [A25]  show  the 
size  of  y  sites. 
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